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ABSTRACT 

Cyclic  voltjmograas  ct  alcrodlsk  alaccrodcs  that  account  for  the  affects 
of  solution  resistivity  double-layer  capacitance,  C^,  .  and  electrochemical 
kinetics  have  been  simulated  by  using  the  explicit  Hopscotch  algorithm  and 
a  conformally  mapped  space-grid.  The  latter  procedure.  In  particular, 
provides  an  efficient  means  of  simulating  voltammograms  at  slower  potential 
scan  rates  where  radial,  rather  than  linear,  dlffiislon  predominates. 
Simulated  voltammetrlc  waveshapes  as  a  function  of  scan  rate  are  reported 
for  systematic  variations  In  pq,  C^j^,  and  the  standard  electrochemical  rate 
constant  k,  ,  both  separately  and  together.  The  waveshape  description  utilizes 
both  the  anodic-cathodlc  potential  peak  separation,  AE^ ,  and  the  halfwave 
potential  for  the  forward  scan  in  relation  to  the  formal  potential, 
(E^-  Ef).  The  latter  waveshape  parameter  Is  suitable  for  analyzing 
voltammograms  where  the  Influence  of  radial  diffusion  Is  sufficient  to  yield 
largely  sigmoidal,  rather  chan  peaked,  waveshapes.  Vhlle  the  use  of  high 
scan  rates  enhances  the  effect  of  electrode  kinetics  upon  the  voltammetry, 
the  deleterious  coupled  Influence  of  fiQ  and  is  also  maximized  under  these 
conditions,  especially  for  larger  values.  At  slower  scan  rates,  however, 
rsslsclve/capacltanee  distortions  are  often  attenuated  greatly,  while 
substantial  kinetic  Information  remains  due  to  the  enhanced  mass  transport 
assocleted  with  radial  diffusion.  The  examination  of  voltamaetrlc  waveshapes 
over  a  wide  range  of  scan  rates,  together  with  proper  consideration  of 
resistance/capacitance  effects,  therefore  constitutes  a  desirable  tactic 
for  evaluating  rapid  electrode  kinetics  using  microdisk  electrodes. 
Representative  simulation  results  are  provided  as  a  guide  to  such  procedures, 
and  as  a  means  of  assessing  Che  degree  of  error  Involved  If 
reslstance/capacltance  effects  ere  neglected. 


Appl^cationa  of  cyclic  vcltamneCry  at  alcrodlsk  alectrodcs  for  the  evalu* 
atlon  ox  the  rates  of  heterogeneous  charge -transfer  reactions, as  veil 
as  of  homogeneous  chemical  processes  have  become  common  in  recent 
years.  The  minute  double-layer  capacitance  seen  at  microelectrodes 
facilitates  the  use  of  extremely  high  potential  scan  races  without  incurring 
extreme  distortions  of  the  voltammogram.^'^*^^  In  the  slow  scan  limit, 
a  well-defined  steady  state  response  can  be  seen,  whereby  the 
diffusion  layer  is  s\ibscancially  larger  chan  the  electrode  dlmenslotis.  This 


yields  relatively  efficient  hemispherical  diffusion  to  and  from  the  electrode 


surface . 


10,11,21-28 


Perhaps  the  most  widely  heralded  feature  of  microelectrode  voltammetry 

is  Che  smaller  degree  of  distortion  due  Co  solution  resistance  as  compared 

29 

with  convantlonal-sizad  electrodes.  This  property  enables  electrochem¬ 
ical  reactions  Co  be  examined  voltammecrically  in  extremely  resistive 
media. However,  it  should  be  stressed  that  while  the  use  of 


microelectrodes  can  reduce  the  severity  of  resistive  distortions,  they  cannot 
eliminate  them.  It  is  therefore  important  to  examine  carefully  the  remaining 
resistance -induced  distortions  in  microelectrode  voltammetry.  This  is  es¬ 
pecially  important  in  the  measurement  of  electrochemical  kinetics,  where 
Che  distorting  effects  of  solution  resistance  can  easily  be  attributed  mis¬ 
takenly  Co  slow  electrode  kinetics. 

These  effects  of  solution  resistance  have  been  modeled  previously  assum¬ 
ing  semi -infinite  linear  diffusion.^’^®  At  microelectrodes,  however, 
Che  attainment  of  this  diffusion  condition  requires  relatively  high  voltam- 
mecrlc  scan  rates.  Resistive  and  migratory  effects  at  microelectrodes  under 
steady-state  conditions,  acces-ed  at  slow  potential  scan  rates,  have  also 
been  examined. Vhile  semi  analytic  solutions  are  available  for 
these  limiting  regimes,  digital  simulation  techniques  are  obligatory  in  the 


int«ra«dlata  c«««.  Halnza  considarad  pravloualy  cha  affacts  of  flnlca 

21  22 

alactroda  kinetics  for  the  latter  circumstances.  ’  Aokl  at  al  have  In- 

23 

vestlgated  linear  sweep  voltammetry  at  microdisks  In  the  reversible  case, 

and  Zoskl  and  co-vorkars  have  developed  new  analysis  methcds  of  "near-steady* 

43 

state"  voltammograsu .  However,  the  effects  of  solution  resistance  and  elec¬ 
trode  capacitance  on  voltammetrlc  waveshapes  In  the  Intermediate  region  have 
not  been  addressed. 

44 

In  a  recent  prallmln^try  note,  we  employed  the  explicit  Hopscotch 

45  46  47-49 

algorithm  ’  as  developed  for  electrochemical  applications  to  simulate 

cyclic  voltammetrlc  waveshapes  at  microdisk  electrodes.  Ue  examined  the 

combined  effects  of  resistance  and  electrode  kinetics  upon  voltammetrlc  peak 

separations.  In  the  present  article,  we  utilize  the  "conformal  map"  algorithm 

50 

for  disk  electrodes  developed  by  Michael  et  al  to  present  a  more  compre¬ 
hensive  and  detailed  report  on  this  topic.  Our  major  objective  Is  to  provide 
a  means  of  diagnosing  and  assessing  quantitatively  the  effects  of  solution 
resistance  and  electrode  capacitance  upon  the  evaluation  of  standard  electro¬ 
chemical  rate  constants,  k^ ,  using  microelectrode  cyclic  voltammetry. 

Computational  Methods 

The  results  presented  In  this  paper  were  genereted  by  a  finite  difference 
algorithm  modeling  the  diffusion  to  a  disk  electrode  Inlaid  flush  with  a 
planar  insulator  during  a  cyclic  voltammetrlc  experiment.  Thus,  the  diffusion 
equation  that  must  be  solved  is 

1  dc  a*Cj  1  dc^  a*Cj 

dt  ar^  r  ar  az* 

whore  r  and  z  are  the  directions  perellel  to  and  norxul  to  the  alect'’ode 
surface,  and  c^  and  are  the  concentration  and  diffusion  coefficient  of 
redox  species  i.  All  results  assume  semi- infinite  diffusional  fields,  with 
Che  reference  electrode  placed  at  infinity.  The  electron- transfer  process: 
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Ox  +  iM  »  Rad  (2) 
la  aodalad  at  the  electrode  eurfece.  The  diffusion  coefficients  of  Ox  end 
Red  ere  essuaed  here  to  be  equal.  The  initial  concentration  of  Ox  is  zero, 
the  solution  containing  only  Red;  with  the  initial  potential  seen  being  in 
the  positive  direction. 

The  slauletlon  can  be  iapleoantad  in  cylindrical  coordinates  with  ex¬ 
ponentially  expending  space  grids  along  each  axis.  This  involves  solving 
siaultaneously  the  flux  aquation  for  the  first  layer  of  solution  above  the 
electrode  surface  and  the  Butler-Volaer  kinetic  relation  (Eqs.  3  and  4, 
respectively) : 

-  -  (CJ-^  .  c;-“)/Az  (3) 

az 


-  Xq  exp<.a  g  AE)  -  exp  ((l-o)g  EA)] 


where  o  is  the  cathodic  transfer  coefficient  (taken  here  as  0.5),  AE  is  the 
difference  between  the  applied  and  the  fomal  potential,  and  Xg  is  propor¬ 
tional  to  tlie  standard  hetergeneous  rats  constant.  Vhen  solved  simultaneously 
for  the  flux,  we  obtain 

«  .  g  ae)  c;-"  -  «»pKi-.)  g  c;-»i  ^ 

”  x;*  ♦  4Z  <«xp(.a  §  4E)  -  g  IE)) 

This  relation  is  than  used  at  each  time  step  to  determine  the  flux  to  the 
electrode  surface.  (For  mors  details  see  ref.  47).  The  two  dimensional 
nature  of  this  problem  requires  a  very  efficient  algorithm  to  avoid  expending 
extrema  amounts  of  computer  time.  The  algorithm  used  initially  here  was  the 
explicit  'fast'  Hopscotch  method. 

Such  direct  solution  of  Eq.  (1)  even  by  the  computationally  efficient 


HopscotcL  algorithm  unfortunately  becomes  cumbersome  when  nonlinear  diffusion 
predominates .  However ,  Kiehael  at  al^^  have  shown  that  a  coordinate  transform 


4 


to  a  spaco  conforalng  to  tha  ahao*  of  tha  diffusion  laysr  at  tho  alactroda 
("conforaal  nap")  and  subsaquant  solution  by  Hopscotch  results  in  ar  algorithm 
which  can  ba  draaatlcally  aora  afficiant  chan  avalxiatlng  Eq.  (1)  in  cylindri¬ 
cal  spaca.  Most  data  presented  here  were  obtained  by  naans  of  the  conformal 
map  procedure,  although  Che  conventional  Hopscotch  method  outlined  above 
was  found  to  be  valuable  for  checking  its  reliability.  The  former  procadoce 
was  iaplaaanted  as  follows.  Equations  (6)  and  (7)  are  the  rsls<’iDns  for 
transformation  from  cylindrical  coordinates  to  tha  conformal  space,  end  Eq. 
(6)  is  Che  conformal  transform  of  Eq.  (1): 

r  -  cos  (#)  cosh  (T)  (6) 

z  -  sin  (#)  sinh  (F)  (7) 

tan  (9)  ^  +  canh  (D  ^  (8) 

ar . 

and  <f  -  (Dj  RT/nF»/a*)**  (9) 

where  a  is  the  radius  of  tha  alactroda,  w  is  Che  scan  rate,  and  9  and  T  are 
tha  axes  of  the  confortel  space. Our  implemencation  of  this  simulation 
uses  the  sane  dimensionless  variables  and  boundary  conditions  described  in 
ret. 50  with  the  addition  of  parameters  for  solution  resistance  and  double- 
layer  capacitance. 

Tha  dimensionless  parameters  describing  solution  resistance,  R, ,  and 
double-layer  capacitance.  are  given  in  Eqs.  (10)  and  (11),  reapsctively: 


dc.  -  <Tj*(ain*«  +  slnh*r]*^  ^  3^c, 

dt  ar* 


/»'  -  4nFaDc*  (nF/RT)  R,  (10) 

7'  -  C,i/(4nFaDc*)  (11) 


These  quantities  are  ralatad  to  the  mors  covon  dimensionless  parametsrs 
38  39 

p,  y  '  employed  in  linear  diffusion  probleau  by 

p  -  p'w/iUc)  (12) 


7  -  Ay'o/a 


(13) 


In  ch«  court*  of  our  ut*  of  chit  Approach,  w*  hav*  nada  cvo  uteful 
obtarvaciont  about  cha  nathod  ictelf.  In  raf.  50.  it  vat  nocad  chat  tha 
■ax  1  arxjra  r-valua,  .  nacattary  to  contain  cha  diffusion  layer  viuhin  the 

tiaulacional  grid  was  detarmined  empirically.  We  have  found  that  this  need 
not  be  Che  case.  The  appropriac*  value  of  can  be  determined  readily  from 

r.„  -  In  (X  -*•  (x»  -  1)‘‘)  (U) 

■her*  X  -  1.0  +  6.0  (15) 

and  where  r  it  the  dimantlonleat  time  of  the  experiment.  These  relations  can 
be  understood  by  noting  that  Eq.  (15)  it  a  modified  version  of  the  traditional 
•6.0(Dt)**  rule*  which  it  used  to  determine  the  thickness  of  the  diffusion 
layer.  Tha  first  term  on  Che  rlghc*hand-tida  of  £q.  (15)  accounts  for  the 
electrode  radius,  which  has  a  value  of  unity  in  dimensionless  cylindrical 
coordlnatet.  Equation  (14)  then  it  simply  a  rearrangement  of  Eq.  (6)  where 
#  -  0,  l.e.  .  at  the  edge  of  the  electrode.  We  have  found  that  when  using 
values  generated  by  Eqs.  (14)  and  (15),  the  simulation  generally  utilizes 
95*98%  of  its  allotted  space. 

The  second  observation  was  made  while  checking  Che  accuracy  of  the  al* 

gorlchm.  To  do  this,  we  input  a  small  a  value,  1  x  10'*,  so  to  approximate 

semi- inf inice  linear  diffusion.  The  result  obtained  with  a  grid  size 

of  20  X  200  eloments ,  as  suggested  in  raf.  50,  did  not  agree  with  that  obtained 

51  52 

froa  coTTvantlonal  theory.  ’  We  found  that  as  o  **  0,  tha  resolution  along 
Cha  r-axis  must  be  increased  in  order  to  faithfully  represent  the  diffusion 
profile.  Figure  3  of  raf. 50  serves  to  explain  further  this  point.  In  fact, 
for  ff  •  1  X  10'*,  the  grid  siza  necessary  to  observe  a  result  within  1%  of 
Chat  obtained  assuming  linear  diffusion  was  20  x  4000.  Clearly,  when  simu¬ 
lating  volCammograms  with  small  o-values,  car*  must  be  taken  to  Insure  cor¬ 
rectness  of  the  result,  and  definitely,  for  very  small  values  of  a,  one  would 
be  foolish  to  us*  this  method.  In  this  regime,  however,  tha  more  conventional 
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•  illation  aachoda  ara  raadlly  appllcabla. 

Siaulatlon  oi  volcamsograas  whan  solution  raslstanca,  ,  is  non-zaro 
is  soaawhat  aore  complex.  Under  these  conditions,  the  passage  of  current, 
i,  generates  a  variable  potential  drop,  IK, ,  within  the  solution  so  that 
the  potential  across  the  double  layer,  no  longer  varies  linearly  with 
tiaa.  The  consequences  of  solution  raslstanca  ara  significant  froa  a 
siaulatlonal  standpoint;  for  now  both  potential  and  currant  are  unknown  and 
intardapandant .  Thus,  a  root-finding  algorithm  visr  be  employed  to  determine 
the  value  of  consistent  with  the  currant  being  passed.  In  our  modifi¬ 
cation  of  the  conformal  sap  to  account  for  resistance,  several  different 
algorithms  were  tested.  A  simple  binary  bisector  worked  consistently ,  but 
increased  computational  time  unacceptably.  The  fastest  root-finding  algorithm 
tested,  a  Newton -Raphs on  procedura^^,  suffered  from  stability  problems  in 
the  presence  of  large  resistances  and  capacitances.  The  Van  Vljngaarden- 
Dekker-Brent  algorithm^^,  however,  provided  sufficient  stabilicy  and 
reasonable  computational  efficiency,  and  was  used  in  all  simulations  presented 
hare  that  involve  non-zero  resistance. 

Three  different  computer  systems  were  used  to  generate  the  results  in 
this  study.  A  Zenith  603^6-based  computer  employing  NDP  Fortran-386,  a 
Digital  Equipment  Corporation  VAX  8800  running  Ultrlx,  and  a  COC  Cyber- 205 
utilizing  VSOS  were  all  employed  hare.  All  calculations  vers  performed  using 
64-blc  precision.  Identical  test  programs  vara  processed  on  all  three  systems 
to  check  machine -independence  of  the  results. 

RESULTS  AHD  DISCUSSION 

ElBgKttdt.KlMUai 

The  most  convenient  and  common  means  of  evaluating  rata  constants  for 
qiiasi-reverslble  electrode  reactions  using  cyclic  voltammetry  involves 
measuring  the  cathodic-anodic  peak  separations,  ,  as  a  function  of  the 
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potential  scan  rate,  ’  Under  the  condition  of  semi -infinite  linear 

diffusion.  for  inoseasurably  large  values  will  be  Independent  of  v 

and  equal  to  2.2  RT/nF.^^  The  observation  of  A£^  values  that  are  larger 

and  that  increase  with  increasing  scan  rata  provides  a  simple  and  well-known 

means  of  evaluating  at  least  in  the  absence  of  significant  distortion 
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due  to  resistive  effects. 

At  a  microdisk  electrode,  however,  the  variation  in  experimental  time¬ 
scale  associated  with  changes  in  sweep  rate  serves  to  modify  the  nature  of 
dlffuslonal  transport  to  the  electrode.  At  slow  scan  rates,  dlffusional 

transport  is  predominantly  radial,  yielding  the  kno%m  steady  state  response 

29 

of  microdisk  electrodes.  At  sufficiently  short  timescales  (i.e.,  at  fast 


scan  rates),  microdisks  behave  identically  to  mors  conventionally  sized  elec¬ 
trodes  exhibiting  linear  diffusion.^  This  varying  dlffusional  character 
at  microdisk  electrodes  causes  the  value  of  A£^  to  depend  significantly  upon 


scan  rata  even  in  the  absence  of  sluggish  electrode  kinetics  or  solution 


resistance . 


21,22,50 


The  situation  in  the  presence  of  finite  electrode  kinetics,  yet  in  the 
absence  of  solution  resistance,  is  illustrated  in  Fig.  1,  in  the  form  of 
plots  of  AE^  versus  log  (•/').  The  results  displayed  were  simulated  assuming 
tha  following  conditions:  a  -  5  x  10'*  cm,  (dashed  traces),  5  x  10'’  cm 
(solid  traces),  D  -  1  x  10'*  cm^  s'* ,  n  -  1,  T  -  298K.  The  sequence  of  curves 


displayed  for  both  electrode  radii  correspond  to  differing  values  of  k, . 
The  lowest  curves  in  each  set  refer  to  the  "reversible  limit*,  i.e.,  for 
k,  -•  •,  and  "he  highest  curves  to  k,  -  0.2  cm  s'*  (see  the  figure  caption 
for  other  details).  At  a  given  sweep  rate,  AE^  increases  as  k,  decreasss- 
a  well  known  result.  Less  familiar,  however,  are  tha  progressiva  increases 
in  AE^  seen  with  decreasing  scan  rats  for  smaller  v  values,  both  in  the  pre¬ 
sence  and  absence  of  finite  electrode  kinetics.  This  tendency  'eflects  the 
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Increasing  efficiency  of  radial  diffusion  as  v  deci  .s. 

Such  scan  rate -dependent  dlffusional  character  is  an  Imporcanc  considera¬ 
tion  for  the  evaluation  of  electrode  kinetics.  One  need  only  glance  at  the 
right-hand  side  of  the  plots  for  the  larger  electrode  (dashed  curves)  in 
Figure  1  to  discern  this.  Here  there  is  a  ca  30  mV  displacement  in  AE^  be¬ 
tween  the  curves  for  the  slowest  rate  constant  (0.2  cm  s'^,  circles)  and 
the  reversible  case  (asterisks).  If  non-linear  diffusion  to  the  electrode 
surface  is  absent  (as  for  the  open  circles  and  squares  in  Fig.  1),  the  cor¬ 
responding  values  of  AE^  differ  by  less  chan  10  mV  at  such  slow  scan  rates. 
Thus,  as  the  values  of  v  decrease  into  the  regime  where  radial  diffusion  is 
Increasingly  prevalent,  the  concurrent  increase  in  dlffusional  efficiency 
enhances  Che  m.<asurabilicy  of  a  given  rate  constant.  Indeed,  this  phenomenon 

is  the  basis  for  the  evaluation  of  electrode  kinetics  with  steady-state  micro- 
10  25 

disk  voltammetry.  '  '  Particularly  noticeable  in  Fig.  1  is  the  markedly 
larger  effect  of  electrode  kinetics  upon  AE^  for  the  smaller  versus  larger 
electrodes  at  slow  scan  races,  arising  from  the  more  efficient  radial  dif¬ 
fusion  in  the  former  case.  At  fast  scan  rates  (i.e.,  the  left-hand  side  of 
Fig.  1),  however,  the  AE^  values  for  Che  two  electrode  sizes  approach  each 
other,  reflecting  the  approach  to  the  linear  diffusion  limit  where  the  effect 
of  electrode  size  on  Che  mass  transport  characteristics  will  disappear. 

Figure  2  displays  the  same  type  of  information  as  Fig.  1,  but  in  a  dimen¬ 
sionless  format  so  to  be  applicable  to  a  wider  variety  of  experimental  situ¬ 
ations.  In  this  figure,  Che  dimensionless  peak  separation  as  given  by 

A<p  -  (nF/RT)  AE^  (16) 

is  plotted  against  log  o,  where  a  is  the  dlffusional  parameter  defined  by 
£q.  (9).  The  dimensionless  quantity  a  provides  a  convenient  yardstick  of 
the  diffusion- layer  thickness  in  relation  to  Che  slecCrode  size.  For  a  «0.1, 
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Ch«  foraar  la  aufflciantly  aaallar  than  cha  lattar  ao  that  aaal-lnfinlca 
llnaar  diffiialon  la  approached.  Incraaalng  a ,  aa  brought  about  by  dacraaalng 
V  and/or  a,  anhancas  tha  Isportanca  of  radial  diffusion  to  tha  overall  aaaa 
transport.  For  example,  at  a  large  electrode,  say  2  ma  radius.  It  would  be 
necessary  to  use  scan  rates  slower  than  ca  0.1  mV  a*^  to  observe  significant 
radial  diffusion  character  (assuming  D  -  10'^  cm  a*^};  clearly  this  Is  pre¬ 
vented  by  convective  mass  transport  at  such  long  timescales.  Decreasing 
tha  electrode  radius  to  2  x  10'*  cm,  however,  will  yield  significant  radial 
diffusion  affects  at  scan  rates  below  100  Vs**,  well  within  the  usual  range 
of  voltammctric  timescales. 

The  bottom  curve  In  Fig.  2  refers  to  the  "revexslble'*  case,  (i.e.  where 
->  •) .  An  empirical  regression  fit  to  this  curve  Is  given  by 


Afp  -  <5.44  ±  0.084)0  +  (2.21  ±  0.027) 


(17) 


21 

These  results  agree  well  with  Heinze's  earlier  simulations.  Note  that 

dfp  Increases  progressively  for  larger  o  values.  For  o  ^  0.5,  however,  tha 

cathodic  and  anodic  voltammetrlc  peaks  are  sufficiently  broad  and  Ill-defined 

aa  to  be  barely  dlscemable.  Eventually  a  "steady-state”  response  Is  at- 

27  28 

talnad,  consisting  of  the  well-known  sigmoid  current -potential  curves.  ’ 

The  nine  upper  curves  In  Fig.  2  were  calculated  for  a  sequence  of  decreasing 

electrochemical  reactivities,  aa  gauged  by  the  conventional  "heterogeneous 
53 

kinetic  parameter” 


(P./D,)*^\ 

[Dir(nF/RT)w]*’ 


'  <18) 


where  0  and  D  are  the  diffusion  coefficients  of  the  oxidized  and  reduced 

O  k' 

redox  forms. 

The  virtue  of  the  dimensionless  representation  in  Fig.  2  Is  that  the 
numerical  results  are  applicable  to  a  range  of  electrode  radii  and  scan  rates. 
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Sine.  ..eh  curv.  corr..pond.  to  .  fix.d  4  v.lu.,  eh.  ratio  of  k,  to  r.«.ln. 
eon.t«nc.  As  th.  dst.  s.t.  in  Fig.  2  sr.  only  .  snsll  r.pr.s.nt.tion  of  .n 
Inflnit.  fwily  of  curv.s,  .  conv.ni.nt  a.thod  of  dotoraining  th. 
v.lu.  of  ^  corr. .ponding  to  given  v.lu..  of  and  of  a  would  b.  of  a.rit. 
W.  h*v.  found  that  ..eh  of  th.  curve,  in  Fig.  2  can  be  fit  to  r.a.onable 
accuracy  by  an  .quation  of  th.  fora  of; 

df,  -  Cj  exp  (C,  logj5  (a))  ♦  c,  (19) 
wh.r.  Cj.  C, ,  C,  are  fitting  paraaet.rs  unique  to  each  value  of  Indeed, 
th.  .olid  trace,  in  Figure  2  w.r.  obtained  in  .uch  a  aanner.  Moreover,  th. 
Cj  coefficient,  each  follow  a  rea.onably  .iapl.  function  of  logj^  (^) ,  .s 
in  Eq.  (20-22): 


where 


a*  ^ 

Cj  -  a^  (1  -  exp  (-(a,  +  +  a^^  +  -  ))) 

•**3 

-  U.87  ±  0.1 
«2  -  1-07  ±  0.03 
■j  -  -1.13  ±  0.06 


and  tdiere  0  -  log^^  (^) . 


a*  -  0.65  ±  0.06 


Cj  -  a,  +  a,  -  a^)*  ♦  .J) 


where 


a,  -  2.127  ±  0.007 
a,  -  0.70  ±  0.05 
Mj  -  0.092  ±  0.12 
a,  -  0.51  t  0.05 


uhera 


C,  -  a,  .xp(.j5  0)  ♦ 

*,  -  1.4  ±  0.1 
•lo  “  ±  0-07 


(20) 


(21) 


(22) 


1,  ' 


-  2.04  ±  0.09 


Given  these  three  relstlons  end  experiaental  values  for  o,  one  can 
search  over  the  Interval  of  experiaentally  reasonable  values  of  ^  until  a 
curve  fitting  the  value  at  the  known  o  value  is  obtained,  rendering 
the  desired  ^  and  hence  the  k^  value.  Although  only  approximate,  this  is  a 
convenient  approach  for  estimating  rate  constants  over  a  vide  range  of  o  from 
peak  separations. 

An  interesting  feat'^re  of  the  curves  in  Fig.  2  is  that  their  vertical 
displacement  (Af^)  from  the  reversible  limit  Increases  surkadly  as  a 
increases.  This  reflects  the  increased  efficiency  of  mass  transport,  and 
hence  a  proportionately  greater  influence  of  electrochemical  kinetics  upon 
the  voltammetrlc  response,  in  the  increasing  presence  of  radial  diffusion. 
Similar  findings  have  been  obtained  by  other  workers. The  results 
in  Fig.  2  therefore  suggest  that  the  most  advantageous  method  of  evaluating 
electrochemical  kinetics  with  microelectroda  cyclic  voltammetTy  would  be  to 
choose  conditions  so  to  maximize  o,  thereby  making  the  diffuslonal  transport 
faster  relative  to  the  electrode  kinetics.  It  is,  alas,  not  that  simple. 
As  already  noted,  the  "peakedness”  of  the  voltasaaogram  also  decreases  under 
these  conditions,  making  it  increasingly  difficult  to  identify  accurately 
the  peak  potentials  from  experimental  data.  This  point  is  illustrated  in  Fig. 
3,  which  consists  of  simulated  voltammograms  for  several  values  of  o,  as 
indicated. 

As  an  sltamatlva  to  peak  potential  separations,  the  half-wave  potential 
E),  for  the  forward  scan  (i.e.,  the  potential  at  which  the  current  reachas 
one  half  the  maximum  value)  can  be  utilized.  Vhile  the  separation  between 
the  E^  values  for  the  forward  and  reverse  scans,  AE{,,  nay  be  employed  in  a 
fashion  similar  to  AE^  values,  one  limitation  is  that  A£)j  approaches  zero 
in  the  steady-state  limit.  A  simple  alternative  is  to  use  Itself,  most 
conveniently  for  the  forward  potential  sweep,  in  relation  to  the  formal  poten- 


tl«l  Ej .  Th«  dlffsranca  (£},•  E^)  c«n  b«  utillz«d  ac  a  klnatlcally  aanaltiva 
paraaatar  chroughoue  tha  regloa  whara  radial  diffusion  la  Important  Including 
tha  staady-stata  limit. 

Figura  4  contains  plots  of  (E^-E^)  varsus  log  (t/'**)  for  tha  aama  set  of 
values  as  In  Fig.  1.  The  x-axia  is  extended,  however,  to  smaller  scan 
rates  [l.a.,  larger  log  values]  so  to  yield  conditions  approaching  the 

steady-state  limit.  Inspection  of  Fig.  4  shows  that  the  sensitivity  of  the 
voltammetry  to  electrode  kinetics,  as  discerned  by  tha  displacement  of  (E^- 
Ef)  for  a  finite  k^  value  above  the  reversible  (k^  -*  ••)  curve,  survives  even 
to  small  scan  rates.  Similarly  to  Fig.  1,  this  results  from  the  offsetting 
Influences  of  the  longer  timescales  and  more  efficient  mass  transport  associa¬ 
ted  with  smaller  scan  rates.  Figure  5  displays  tha  same  type  of  Information, 
but  In  a  dimensionless  format,  as  a  plot  of  (Ctj*  (f)  versus  log  o,  where 
(V  (f)  -  (nr/RT)(Ei,-E,)  (cf  Fig.  2). 

Sgiutim  Rtalatflpfig 

The  presence  of  solution  resistance  has  long  been  known  to  cause  signifi¬ 
cant  distortions  to  cyclic  voltammo grams,  most  notably  to  enlarge  tha 
1  38  39 

values.  ’  ’  Unfortunately,  the  functional  nature  of  this  influence  upon 

is  often  similar  to  that  brought  about  by  finite  electrode  kinetics.  The 

unambiguous  Identification  of  kinetic  and  resistive  effects  In  their 

simultaneous  presence  Is  therefore  clearly  lisportant.  While  this  Issue  has 

1  3 

been  the  focus  of  several  studies,  ’  these  have  dealt  only  with  the  case 

of  pure  linear  diffusion.  Below,  we  provide  some  examples  of  the  offsets  of 

resistance  on  voltammogrsms  also  Influenced  by  the  changing  dlffuslonal 

character  of  microdisk  electrodes. 

54 

Newman  has  shown  that  the  resistance,  ,  for  a  disk  electrode  in 
isotropic  media  with  tha  auxiliary  and  reference  electrodes  an  Infinite 
distance  away  is  given  by 
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R,  -  Pfl/A*  (23) 

whars  Pq  is  the  solution  resistivity.  Although  Eq.  (23)  shows  that  R,  will 
increese  es  the  electrode  redius  decreeses,  this  effect  is  oore  then  offset 
by  the  diminished  currents  obtained  under  these  conditions,  yielding  smeller 
IR^  potential  drops.  The  currant  under  conditions  of  linear  diffusion  (a  - 
0)  is  proportional  to  the  electrode  area  (and  hence  a*),  so  that  IR,  will  de¬ 
crease  linearly  with  decreasing  a.  At  the  extreme  of  steady -state  diffusion 
23 

to  a  microdisk,  (a  -*  «•) ,  the  limiting  current  becomes  proportional  to 

the  electrode  radius,  whereupon  the  IR,  drop  is  independent  of  electrode 
41  42 

area  (vide  infra).  '  At  finite  non-zero  values  of  o,  an  intermediate 

situation  will  prevail. 

While  it  has  been  common  to  n£glect  resistive  effects  in  kinetic  applica¬ 
tions  of  microelectrode  voltammetry,  this  is  clearly  a  dangerous  tactic. 
This  point  is  Illustrated  straightforwardly  in  Fig.  6,  as  plots  of  AE^  versus 
log  (u  ^)  for  various  p^  values  but  in  the  absence  of  finite  electrode 
kinetics  (i.e.,  for  k,  ■*  •) .  As  before,  the  solid  and  dashed  traces  refer 
to  0.5  and  5  pm  radius  electrodes,  respectively.  The  values  ere  selected 
so  to  span  the  values  commonly  encountered  in  organic  solvent  systems.  Thus 
the  increasing  pq  values  in  Fig.  6  (from  200  to  2670  0  cm)  correspond  to 

nltromethane ,  formamide,  methylene  chloride,  ethanol,  and  tetrahydrofuran, 
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each  containing  0.1  tetrabutylammoniua  perchlorate  (see  figure  caption 
for  further  details).  The  lowest  curve  in  each  sat  corresponds  to  p^  -  0. 

Comparison  of  Fig.  6  with  Fig.  1  shows  that  the  effect  of  solution  resis¬ 
tance  upon  AE^  as  a  function  of  scan  rate  is  qualitatively  similar  to  that 
observed  for  finite  electrode  kinetics.  Nevertheless,  there  are  two  note¬ 
worthy  differences.  Most  strikingly,  the  effect  of  altering  the  electrode 
size  at  a  given  scan  rata  is  qualitatively  different.  While  the  influence 
of  solution  resistance  upon  AE^  diminishes  as  the  electrode  radius  decreases 


(Fig.  6),  tha  rola  axartad  by  alactroda  klnatlca  la  anhancad  under  chaae 
condlciona  (Fig.  1).  Aa  noted  above,  the  latter  a  'act  becomea  Increaalngly 
latportant  aa  the  electrode  alze  decreaaes.  and  hence  the  degree  of  radial 
diffusion  Increases.  In  addition,  careful  comparison  of  Figs.  1  and  6  reveals 
that  even  at  the  larger  electrode  the  Influence  of  electrode  kinetics  upon 
i.ecoaes  greater  relative  to  the  corresponding  effect  of  solution  resis¬ 
tance  aa  the  scan  rate  decreases.  This  phenomenon  also  has  its  origin  in 
the  enhanced  mass  transport  efficiency  that  ensues  as  the  degree  of  radial 
diffusion  is  increased. 

These  comparisons  also  serve  to  eriphaslze  the  central  role  played  by 
the  electrode  size  in  Influencing  the  measurability  of  electrochemical  kin¬ 
etics  in  the  face  of  distortions  from  solution  resistance.  The  common  tactic 
of  increasing  the  scan  rate  can  enhance  greatly  the  effect  of  electrode  kine¬ 
tics  upon  AEp ,  particularly  if  the  degree  of  radial  diffusion  is  small  (i.e., 
the  left-hand-side  of  Fig.  1).  However,  the  effect  of  solution  resistance 
upon  also  can  Increase  markedly  under  these  conditions ,  tending  to  negate 
the  enhanced  effect  of  electrode  kinetics.  By  examining  A£^  at  lower  scan 
rates,  on  the  other  hand,  the  effect  of  electrode  kinetics  relative  to  solu¬ 
tion  resistance  can  be  enhanced  dramatically  if  sufficiently  small  electrodes 
are  employed  so  to  provide  efficient  radial  diffusion. 

Similar  arguments  apply  if  half-wave  potentials  rather  than  peak  separa¬ 
tions  are  employed.  Figure  7  is  a  plot  of  (E^-  E, )  versus  log  (w*^)  for  the 
same  set  of  solution  resistance  values  and  other  parameters  as  in  Fig.  6. 
The  effect  of  resistance  is  again  seen  to  be  minimized  at  low  scan  rates  as 
wall  as  fer  smaller  electrodes.  One  additional  point  is  evident  tjpon  com¬ 
paring  Fig.  7  with  the  corresponding  plot  for  electrode  kinetics  in  Fig.  4. 
Vhlle  the  sensitivity  to  electrode  kinetics  is  maintained  over  a  vide  range 
of  scan  rates  in  the  radial  diffusion  regime,  the  effect  of  resistance  (if 
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■Ignlflcsnt)  diiplays  •  similar  Invarianca.  Consaquently,  than,  one  must 
ba  careful  to  estimate  the  degree  to  which  the  known  solution  resistance 
should  distort  experimental  voltammcgra/is  before  asserting  chat  and  E^ 
values  that  deviate  from  the  "Ideal  re/ersible”  limits  signal  the  presence 
of  measurable  electrode  kinetics. 

With  the  aim  of  providing  a  more  complete  numerical  picture  of  Che  ef¬ 
fects  of  solution  resistance  on  microdisk  cyclic  voltammetry,  Figs.  8  and  9 
display  plots  of  and  versus  log  o  for  a  series  of  values  of 
the  dimensionless  resistance  parameter,  p  [Eqs.  (10),  (12)].  A  cursory 
examination  of  both  figures  reveals  that  the  separation  between  the  uppermost 
curve  (p  •  10)  and  Che  lowest  curve  (p  •  0)  increases  as  a  increases,  l.e., 
from  left  to  right.  This  Is  caused  by  Che  relative  increase  in  current  magnl* 
tuda  as  a  increases,  l.e.,  by  Che  more  efficient  diffusion  observed  near 
steady- state  conditions.  The  dimensionless  format  la  deslgnod  to  maintain 
Che  resistive  drop  constant  relative  to  the  scan  rata  in  the  linear  diffusion 
limit.  As  a  result  of  the  Increased  dlffuslonal  efficiency  encountered  at 
lover  scan  rates,  Che  magnitude  of  the  resistive  drop  Increases  as  the  steady- 
stats  limit  Is  approached.  Fortunately,  one  Is  unlikely  to  encounter  a 
situation  where  p  Is  large  for  steady-state  conditions,  except  for  extremely 
high  solution  resistances  such  as  are  encountered  In  the  absence  of  supporting 
alactrolyta.  In  this  case,  however,  Oldham  has  pointed  out  that  Newman's 

Isotropic  condition  la  no  longer  valid,  and  the  ohmic  trsatmant  of  resistance 

34 

Is  likavlse  Inappropriate. 

Tha  results  discussed  so  far  neglect  Che  additional  Influence  of  double - 
layer  capacitance,  on  the  voltammetrlc  response.  In  the  absence  of 
significant  effects  from  solution  resistance,  a  constant  non-zero  capacitance 
will  provide  merely  a  background  nonfaradalc  current  without  altering  Che 
voltammetrlc  waveshape.  In  tha  presence  of  solution  resistance  effects. 
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hovavar,  tha  alactroda  capacltanca  can  ylald  additional  voltanunetrlc 
dlatortlona  alnca  the  tlaa  derivative  of  tha  double -layer  potential,  dZ^^/dt, 
and  hence  the  nonfaradalc  current,  can  vary  aubstantlally  with  potential 
even  vhen  is  potential- independent. This  situation  yields  a  close 
coupling  between  the  voltammetric  distortions  arising  from  R,  and  in 

that  tha  presence  of  tha  latter  tends  to  enhance  the  affect  of  the  former. 

This  point  is  illustrated  in  Figs.  10  and  11.  in  the  form  of  plots  of 
and  <Z^-  Ej),  respectively,  versus  log  (v***)  (cf  Figs.  6.7)  for  an  elec¬ 
trode  radius  of  S  tm.  The  open  and  filled  symbols  refer  to  pq  values  of 
725  and  1800  0  cm'^,  respectively.  (These  values  were  chosen  since  they 
are  appropriate  for  voltammetry  in  the  moderately  resistive  media  methylene 

chloride  and  ethanol,  respectively,  both  containing  0.1  H  supporting 
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electrolyte  ).  The  trio  of  points  in  these  two  sets  (circles,  sqvures,  and 
triangles)  refer  to  values  of  5,  15,  and  25  pF  cm*^,  respeerlvely.  These 
plots  show  the  extent  to  which  the  presence  of  larger  values  accentuates 
Che  effect  of  solution  resistance.  At  the  highest  scan  rats  (10*  V  s'*) 
shown  in  Fig.  10,  for  example,  increasing  fivefold,  from  5  to  25  ^F 
cm'^,  increases  to  an  extent  similar  to  enhancing  R,  by  ca  2.5-fold,  from 
725  Co  1800  n  cm'*.  The  relative  importance  of  such  effects,  however, 
diminishes  sharply  towards  smaller  scan  rates  (Figs.  10,  11}  as  expected 
since  the  magnitude  of  the  nonfaradalc  currents  is  decreased  under  these 
conditions. 

These  results  indicate  chat  the  use  of  hi^  scan  rates,  often  sought 
after  as  a  means  of  evaluating  rapid  electrode  kinetics,  can  be  problematic 
for  this  purpose.  Under  conditions  where  solution  resistance  is  known  (or 
suspected)  to  distort  significantly  tha  voltammograms ,  It  will  often  be  neces¬ 
sary  to  account  additionally  for  capacitance  affects.  A  difficulty  with 
the  latter  is  chat  is  inherently  a  surface -sensitive  property,  its  mag- 


nltud*  being  determined  by  the  Interfeclel  structure  composition  (end  even 
by  electrode  construction  techniques'^).  In  addition,  significant  capacitance 
contributions  can  arise  from  electrode  leads  and  connections  ("stray  capaci¬ 
tance").  T.t  Is  therefore  Important  to  evaluate  carefully  the  occurrence 
and  magnitude  of  such  coupled  resistance-capacitance  effects. 

Combination  af. Solution  Resistance.  Double-Laver  Capacitance  and  Electrode 
Kinetics 

Ideally,  one  would  prefer  to  present  a  complete  picture  of  the  effects 
of  resistance,  capacitance  and  kinetics  when  they  are  simultaneously  present. 
Tais  undertaking,  ever.  In  the  linear  diffusion  limit,  would  be  a  formidable 
task.  By  Including  the  scan  rate -dependent  dlffxislonal  character  obtained 
at  a  microdisk.  It  becomes  completely  Impractical!  Ve  nonetheless  wish  to 
provide  a  brief  overview  of  the  coupled  effects  that  all  these  factors  can 
exert  on  vcltammetrle  wave  shapes.  To  that  end,  we  provide  some  specific 
examples  of  the  combined  effects  of  resistance,  capacitance,  and  electrode 
kinetics.  Included  are  some  examples  of  the  errors  that  can  accrue  If  resis¬ 
tive  and  capacitive  distortions  are  Ignored  In  the  evaluation  of  electrode 
kinetics. 

The  first  example  Is  contained  In  Figs.  12  and  13.  These  figures  contain 
corresponding  plots  of  and  (E^  •  E^).  respectively,  against  log  (*'*‘‘). 
hM  before,  the  solid  and  dashed  lines  refer  to  0.5  and  5  (a  radius  electrodes, 
respectively.  Results  are  shown  In  Figs.  12  and  13  for  a  pair  of  solution 
resistances:  100  Q  cm  (^pen  symbols)  and  1000  Q  cm  (filled  symbols).  For 
both  of  these  pq  values,  data  are  shown  for  a  relatively  low  and  high  k^ 
value,  0.2  cm  s*^  (circles)  and  5  cm  s*^  (squares).  The  double-layer  capaci¬ 
tance  is  fixed  in  each  case  at  5  pF  cm*’;  ocher  parameters  are  as  utilized 
above. 

As  expected,  the  voltammetrlc  responses  at  the  larger  electrode  are 


aff«ct«d  xubatanclally  mora  by  solution  rasistsnce  than  chose  for  the  smaller 
electrode,  especially  at  high  scan  rates.  The  degree  to  which  solution  re¬ 
sistance  alters  both  A£^  and  E^,  however,  is  dependent  significantly  on  the 
extent  to  which  electrode  kinetics  also  affects  these  parameters.  Thus  the 
increase  in  obtained  for  a  given  value  by  increasing  pq  (discerned 
by  comparing  corraapcnding  open  and  filled  symbols  in  Fig.  12)  is  noticeably 
lass  for  k,  <•  0.2  cm  s~^  chan  for  5  cm  s' In  other  words,  Che  extent  of 
additional  resistive  distortion  upon  tends  Co  diminish  somewhat  for 

voltammograms  Chat  also  feature  Che  effects  of  more  sluggish  electrode  kine¬ 
tics.  This  coupling  between  the  effects  of  solution  resistance  and  electrode 
kinetics  arises  primarily  from  the  diminution  in  the  currents  Chat  result 
from  slower  kinetics.  Ve  also  pointed  out  this  non-additive  nature  of  solu¬ 
tion  resistance  and  electrode  kinetic  effects  upon  dE^  in  our  preliminary 
44 

communication.  Similar  effects  are  also  seen  when  (E^-  )  is  utilized 

as  Che  voltammetric  parameter  (Fig.  13). 

A  necessary  consequence  of  this  coupling  >'etwaen  resistance  and  kinetic 
effects  is  Chat  Che  additional  effect  of  electrode  kinetics  upon  dE^  and 
(E)j-  Ej)  tends  to  diminish  somewhat  in  the  presence  of  distortions  from  solu¬ 
tion  resistance.  This  situation  emphasizes  the  need  to  select  conditions, 
such  as  slower  scan  rates  and  smaller  electrode  radii,  so  that  the  effect 
of  solution  resistance  is  minimized  or  at  least  characterized. 

One  other  feature  of  Fig.  13  worth  mentioning  here  is  that  the  displace¬ 
ment  in  (Z^-  Ef)  for  each  pair  of  curves  corresponding  to  the  different 
values  become  not  only  insensitive  to  scan  rata  for  small  w,  but  are  also 
independent  of  the  electrode  radius  under  these  conditions.  Although  this 
phenomenon  does  not  appear  to  be  particularly  useful  in  the  present  context, 

it  is  In  agreement  with  analytical  predictions  that  refer  to  the  steady- state 
41  42 

limit.  ’  As  noted  above,  an  additional  complication  wrought  by  the  occur- 
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ranca  of  significant  voltanaacric  distortions  from  solution  raslstanca  is  that 
cha  axtant  of  this  affact  vill  also  ba  dapandant  on  tha  magnltuda  of  the 
doubla-layar  capacitanca.  Figures  lA  and  15  contain  plots  in  tha  same  format 
as  Figs.  12  and  13,  respectively,  that  illustrate  the  effects  of  altering 
tha  double-layer  capacitance,  from  5  n?  cm'*  (circles)  to  25  >iF  cm'* 

(squares).  Each  of  these  curves  in  Figs.  14  and  15  refer  to  a  rate  constant 
of  0.2  cm  s' ^  .  (Tha  other  symbolism  and  conditions  are  as  in  Figs.  12  and 
13). 

Inspection  of  Figs.  14  and  15  shows  that  variations  in  the  double-layer 
capacitance  as  well  as  the  solution  resistance  can  exert  a  significant  in¬ 
fluence  on  both  and  (£^•  )  in  the  presence  of  finite  electrode  kinetics, 

especially  for  larger  electrodes  and  high  scan  rates.  As  before  (Figs.  12, 
13),  tha  effects  of  and  on  these  parameters  are  strongly  coupled. 

Thus  increasing  fivefold,  from  5  to  25  >*F  cm'*  for  the  larger  p^  value 
1000  Q  cm,  increases  at  high  scan  rates  to  an  extent  comparable  to  the 
distortions  incurred  by  such  a  large  solution  resistivity  alone.  Such  effects 
incurred  by  altering  p^  and/or  at  such  large  resistivities  can  easily 
be  comparable  to  those  incurreo  by  substantial  changes  in  .  At  the  highest 
Pq  value  and  scan  rates  in  Fig.  14,  for  example,  the  effect  upon  of  in¬ 
creasing  from  5  to  25  (tF  cm''  can  be  offset  almost  entirely  by  altering 
by  a  similar  factor,  from  0.2  co  1  cm  s'*.  Consequently,  the  reliable 
extraction  of  k^  under  these  conditions  requires  quantitative  information 
on  both  pq  and  .  Such  coupled  affects  provide  a  ciear  warning  to  the 
experimentalist  attempting  to  evaluate  rapid  electrode  kinetics  by  means  of 
microelectrode  voltammetry.  Tha  use  of  high  scan  rates  and  larger  (ca  5  pm 
radius)  electrodes,  employed  commonly  for  this  purpose,  are  precisely  the 
conditions  for  which  the  coupled  distortions  from  solution  resistance  and 
electrrde  capacitance  effects  are  most  severe. 


A  further  illustration  of  the  delsterlous  influence  of  such  factors 
upon  the  evaluation  of  electrode  kinetics  is  ;iven  in  Fig.  16.  Simulated 
voltaamograms  were  generated  at  a  sequence  of  scan  rates  from  1.0  to  10*  V 
for  values  equal  to  0.1,  0.2  and  0.5  cm  s*^  (triangles,  squares,  and 
circles,  respectively,  in  Fig.  16),  with  a  -  5  ttm,  pq  -  200  Q  cm' \  and 
”  10  >iF  cm'*  (This  solution  resistivity,  a  ‘moderate*  value,  approximates 
that  appropriate  for  methylene  chloride  containing  0.5  ti  tetraalkylammonium 
salt).  Apparent  estlmaCes  of  k^ ,  k, (app)  ware  extracted  from  each  of  these 
voltammograms  by  using  the  procedure  outlined  above,  but  assuming  that  p^  • 
0.  This  procedure  therefore  mimics  the  common  assumption  that  solution  re¬ 
sistance  sffa:ts  can  be  neglected  in  microelectrode  voltaanetry.  The  degree 
of  systematic  error  resulting  from  this  neglect  of  solution  resistance  is 
displayed  in  Fig.  16  in  the  form  of  kj(app)/k^  ratios  plotted  versus  log 
{y'^) .  Inspection  of  Fig.  16  shows  that  while  k^ (app)/kj  ■  1  for  y  <  100  V 
s'*,  this  ratio  diminishes  systematically  towards  higher  scan  rates, 
especially  for  larger  k, .  These  results  affirm  that  the  neglect  of  solution 
resistance  in  the  evaluation  of  k^  becomes  increasingly  serious  at  higher 
scan  races,  k^ (app)  estimates  that  are  progressively  smaller  than  the  true 
values  being  recovered  under  these  conditions. 

A  final,  perhaps  more  visual,  example  of  resistive  effects  uoon  the 
evaluation  of  electrode  kinetics  is  contained  in  Figs.  17A  and  B.  The  former 
displays  two  sets  of  three  cyclic  voltammograms,  simulated  at  a  relatively 
high  scan  rats,  2500  V  s'*,  and  for  the  usual  conditions  considered  here, 
with  a  -  5  >flB,  »  10  mF  cm'*  ,  with  solution  resistivities  of  0  (solid 
traces),  130  (dashed  trace.**),  and  720  0  cm'*  (dotted  traces).  The  latter 
two  sq  values  are  appropriate  for  0.1  H  tetrabutylammonlum  perchlorate  in 
acetonitrile  and  methylene  chloride,  respectively.  The  inner  trio  of  voltam- 
mograms  (i.e.,  with  the  smaller  AE^  values)  were  calculated  for  k^  •,  and 
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ch«  outer  trio  refer  to  k,  ■>  0.2  ca  s~^.  Cxeainetion  of  Fig.  17A  shows  thst 
distinctly  different  voltsaaetrlc  weveshspes  ere  obtained  for  e  given 
velue  es  the  solution  resistance  is  altered  over  this  coomonly  encountered 
range.  Indeed,  the  Increased  A£^  and  related  voltametiric  distortions  attend¬ 
ing  such  alterations  in  the  solution  aedlua  could  easily  be  identified  mis¬ 
takenly  as  reflecting  substantial  solvent- induced  variations  in  the  electron- 
transfer  rate  constant. 

Figure  17B  contains  corresponding  cyclic  voltaaaetrlc  data  as  in  Fig. 
17A,  but  for  a  slower  scan  rata,  25  V  s*^.  Only  two  curves  are  shown,  cor¬ 
responding  again  to  k^  -*  •  (solid  trace)  and  k^  -  0.2  ca  s' ^  (dotted  trace). 
Unlike  the  data  at  the  faster  scan  rate,  (Fig.  17A)  the  voltaaaograas  In 
Fig.  17B  slaulated  for  the  three  values  are  essentially  coincident,  i.e., 
the  effect  of  solution  resistance  is  negligible  under  these  conditions. 
The  influence  of  electrode  kinetics  is  nonetheless  still  significant  at  the 
slower  scan  rate,  as  seen  from  the  noticeable  difference  betwoen  the  two 
voltaaaograas  In  Fig.  17B.  Admittedly,  however,  the  degree  of  kinetic  Infor- 
nation  available  under  these  conditions  is  clearly  smaller  than  that  gleaned 
for  the  fast  scan  rate  in  Fig.  17A,  as  discerned  froa  the  markedly  greater 
increase  in  AE^  resulting  froa  decreasing  k,  in  the  latter  case.  Clearly, 
it  is  iaportant  to  select  aeasureaent  conditions  where  the  desired  degree 
of  kinetic  inforaation  is  Buintalned  while  ainiaizlng  (or  at  least  control¬ 
ling)  the  known  influence  of  coupled  solution  resistance/electrode  capacitance 
effects. 

Concluding  Remarks 

This  study  demonstrates  that  the  reliable  slnnilacion  of  cyclic 
voltaaaograas  at  microdisk  electrodes,  even  in  the  presence  of  substantial 
radial  diffusion  and  with  significant  effects  froa  both  solution  resistance 
and  finite  electrode  kinetics,  is  coaputatlonally  feasible  by  using  a  con- 
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fonaal  spAca  grid.  This  slgorlthn  is  sufficisntly  sfficisnt  that  such  siaula- 
tlons  can  be  effectively  generated  on  high  perforaance  personal  coaputers. 
On  faster  computer  systems,  a  variant  of  this  algorltha  could  also  be  used 
CO  fit  simulated  volcamaograaa  to  experiaental  data  in  an  interactive  fashion, 
as  is  coamon  for  the  case  of  linear  diffusion. 

The  present  results  show  that  the  ability  to  evaluate  rapid  electrode 
kinetics  can  be  iapalred  substantially  at  aicrodisk  electrodes  even  in 
aoderately  resistive  solvents.  Vhile  the  coamon  tactic  of  employing  high 
scan  rates  enhances  the  degree  of  kinetic  information  attainable,  the  dele¬ 
terious  influence  of  solution  resistance  can  be  sufficiently  great  under 
these  condicloiis  so  to  impair  severely  the  extraction  of  reliable  rate  con¬ 
stant  values.  This  difficulty  is  made  more  acute  by  the  coupled  influence 
of  double-layer  capacitance,  which  can  act  to  magnify  substantially  the  dis¬ 
torting  influence  of  solution  resistance  at  high  scan  rates. 

In  the  light  of  these  considerations,  in  order  to  extract  reliable  elec¬ 
trode  kinetic  data  it  is  desirable  to  analyze  aicrodisk  voltamaetric  data 
obtained  for  a  range  of  scan  rates  that  also  encompasses  the  regime  where 
radial  diffusion  contributes  importantly  to  the  mass  transport.  Other  than 
at  the  steady-state  limit  the  analysis  of  this  regime  is,  admittedly,  computa¬ 
tionally  quite  demanding.  Nevertheless,  the  greater  mass -transport  efficiency 
engendered  by  radial  diffusion  acta  to  enhance  the  influence  of  electrode 
kinetics  on  Che  voltammetric  waveshape  relative  to  the  undesired  coupled 
effects  exerted  by  solution  resistance  and  double-layer  capacitance.  As  a 
consequence,  Chen,  the  careful  examination  of  voltammograms  over  an  appropri¬ 
ately  large  range  of  scan  rates,  together  with  the  appropriate  knowledge  of 
Pq  and  ,  may  be  considered  to  be  a  prerequisite  for  the  reliable  extraction 
of  electrochemical  rata  constants. 
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ggJRS.  gAfTIOPS 
£1&.  1 

Plots  of  cyclic  volcammetrlc  p«ak  sepsrstlons,  observed  et  e  S  >ui 
rsdius  (dsshed  traces)  and  a  0.5  ^  radius  (solid  traces)  microdisk,  assuming 
D  -  10'^  cm*  s'*  n  -  1,  T  -  298*K,  versus  logCi/"**)  for  several  different 
heterogeneous  rate  constants;  kj(cm'*):  •»,  asterisks;  5,  inverted  triangles; 
2,  diamonds;  1,  upright  triangles;  0.5,  squares;  0.2,  circles.  Open 
symbols  for  k,  “  •  (squares),  and  k^  -  0.2  cm  s'*  (circles). 

ru.  2 

Plots  of  versus  log(o)  for  the  following  values  of  the  dimensionless 
kinetic  parameter  •,  crossed  inverted  triangles;  10,  crossed  upright 
triangles;  5,  crossed  squares;  2,  crossed  circles;  1.  plus  signs;  0.5,  in* 
verted  open  triangles;  0.2,  open  diamonds;  0.1,  upright  open  triangles; 
0.05,  open  squares;  0.02,  open  circles. 

Fig.  3 

Voltammetric  waveshapes  (for  k^  -*  •)  seen  under  different  diffusions! 
conditions  at  a  microdisk,  for  the  specific  a  values  indicated. 

£U^ 

Plots  of  (El,  *  E,)  versus  log(i''^)  for  several  different  rate  constants. 
Symbols  and  conditions  are  as  in  Fig.  1. 

Fig.  5 

Plots  of  (El,  *  Ef)  versus  log(o)  for  a  series  of  values  of  Symbols 
are  as  in  Fig.  2. 

Elt.  $ 

Plots  of  versus  log(w'^)  for  voltammograms  generated  for  k,  **  * 
but  with  nonzero  solution  resistivities,  (0  -  cm)  as  follows:  Pq  ~ 
asterisks;  200,  Inverted  triangles;  450,  diamonds;  725,  upright  triangles; 
1800,  sjiiares;  2670,  circles.  Reactant  concentration  (C,^}  was  2  x  10**  {1. 
Other  conditions  as  in  Fig.  1. 

fig.  7 

Plots  of  (E|,  •  Ej)  versus  log(i/*  )  for  voltammograms  generated  for 
conditions  as  in  Fig.  6. 

fu^a 

Plots  of  versus  log(o)  for  a  series  of  different  values  of  the 
dlmenalonlass  resistive  parameter,  p,  as  follows:  ^  -  0,  asterisks;  0.1, 
closed  circles;  0.2,  open  squares ;  0.5,  inverted  closed  triangles ;  1.0,  closed 
diamonds;  2.0,  upright  closed  triangles;  5.0,  closed  squares;  10,  open 
circles . 
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Plots  of  (f^  -  fj)  versus  log(o)  for  s  series  of  different  values  of 
p.  Syabols  and  conditions  are  as  In  Fig.  8. 

TiS^lQ 

Plots  of  versus  log<v**’)  for  a  S  pa  radius  electrode  with  conditions 
In  Fig.  Z  for  pq  -  725  Q-ca  (open  syabols)  and  1800  Q>cb  (closed  syabols), 
with  Cjj  -  5  pF  ca**  (circles),  15  pF  ca**  (squares)  and  25  pF  ca'* 
(triangles).  Reactant  concentration  (C^)  was  2  x  10*’  M,  diffusion  coeffi¬ 
cient  (D)  -  1  X  10“  ®  ca^  s'^ 

FU^Il 

Plots  of  (E^j  *  versus  log(w***)  for  conditions  as  In  Fig.  10. 

Fix.  12 

Plots  of  versus  log(v*‘’)  for  conditions  as  follows:  a  -  0.5 

pa  (solid  traces),  5  pa  (dashed  traces):  k,  -  0.2  ca  s'^ ,  (circles),  5.0  ca 

s'^  (squares),  p^^  ••  100  Q-ca  (open  syabols),  1000  0-ca  (closed  s3n>bols).  Ocher 

condition/:  -  5  pF  ca**.  n  -  1,  T  -  D  -  10*’  ca*  s**,  C*  -  2  x 

10*  *M, 

FU..  13 

Plots  of  (E)y  *  F{)  versus  log(w*'*)  for  conditions  as  In  Fig.  12. 

Flg>  14 

Plots  of  versus  log(i/*^)  for  conditions  as  follows:  a  •  0.5  pa 
(solid  traces) ,  5  pa  (dashed  traces  );  -  25  pF  ca'*  (squares),  5  pF  ca** 

(circles),  Pq  -  100  Q-ca  (open  syabols),  1000  Q-ca  (closed  syabols).  Other 
conditions:  k,  -  0.2  ca  s**;  -  2  x  10*’M,  0-10*’  ca*  s**,  n  -  1,  T  - 

298’X. 

Fix.  15 

Plots  of  (Etj  *  E{)  versiis  log(w'^)  for  conditions  as  In  Fig.  14. 

Flf.  16 

Ratio  of  "apparent”  Co  actual  rate  constant,  k,  (•Pp)A,i  versus  log 
(v***)  for  the  following  conditions:  a  -  5  pa,  n  -  1,  T  ••  298*K,  pq  -  200  0- 
ca,  -  10  pF  ca'*,  C"  -2  x  10**M,  D  -  10**  ca*  s**,  k,  ••  0.5  ca  s'*  (cir¬ 
cles),  0.2  Cl.  s'*  (squares),  0.1  ca  s'*  (triangles).  The  k,  (app)  values 
vers  extracted  froa  the  alaulated  volcaaaograas  by  ignoring  solution  resis¬ 
tance  effects  (see  text) . 

Flt..JJ 

Cyclic  voltaaaograsu  alaulated  for  k,  -»  «•  and  0.2  ca  s'*  at  sweep  rates 
of  (a)  2500  V  s'*  and  (b)  25  V  s'*.  Ceneral  conditions  are:  a  -  5  x  10'* 
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cm,  D  "  1  X  10**  c**  ■'*,  Cjj  -  10  nT  cm~* ,  -  2  x  10**  H.  In  (*) ,  solution 

raslsClvltias  (pq)  ara:  0  (solid  tracas),  130  (dsshad  tracas),  and  720  0 
ca  (dottad  traces).  Tha  Inner  and  outer  trio  of  voltamograas  refer  to 
-•  «•,  and  -  0.2  ca  a**,  respectively.  In  (b) ,  solid  and  dotted  traces 
refer  to  k,  -*  •  and  k^  -  0.2  ca  a**,  raapactlvaly  (curves  for  all  resistivi¬ 
ties  above  are  essentially  coincident  -  sea  text) . 
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